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Abstract 



We discuss the effect induced on the neutralino-nucleon cross-section by 
the present uncertainties in the values of the quark masses and of the quark 
scalar densities in the nucleon. We examine the implications of this aspect on 
the determination of the neutralino cosmological properties, as derived from 
measurements of WIMP direct detection. We show that, within current theo- 
retical uncertainties, the DAM A annual modulation data are compatible with 
a neutralino as a major dark matter component, to an extent which is even 
larger than the one previously derived. We also comment on implications of 
the mentioned uncertainties for experiments of indirect dark matter detection. 
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I. INTRODUCTION 



The sensitivities of the experiments of direct search for WIMPs have remarkably im- 
proved in recent years, allowing now the exploration of sizeable regions of the physical 
parameter space of specific particle candidates for dark matter. This is the case of the 
neutralino [|TJ, for which some direct detection experiments are already capable of investi- 
gating significant features in domains of the parameter space which are also under current 
exploration at LEP2. 

The signal searched for in experiments of WIMP direct detection is a convolution of the 
WIMP velocity distribution in the halo with the quantity p x <J e u where p x is the local WIMP 
matter density M and a e i is the WIMP-nucleus elastic cross-section. Under the assumption 
that the WIMP has the two following properties: i) its cross-section with matter is dominated 
by coherent effects, ii) its (spin-independent) couplings are essentially the same for protons as 
for neutrons, then one can straightforwardly factor out a scalar (spin-independent) WIMP- 
nucleon cross-section, <J^i^° n \ in a e i. In this instance, the information derivable from any 
experiment of WIMP direct search may be directly formulated in terms of the quantity 
Px^scaiar 6011 '' j an d the sensitivities of various experiments may be easily compared among each 
other ||. Properties i) and ii) are actually satisfied almost everywhere in the supersymmetric 
parameter space when the WIMP candidate is a neutralino, which is the case explicitly 
discussed in the present paper. 

Some of the WIMP direct search experiments are already sensitive, or are on the verge 
of becoming sensitive, to the range p^ 3 ^^ ^ ~ a few • 1CT 9 4- 1 • 1CT 8 nbarn, where p° 3 
denotes the local density normalized to the standard value of 0.3 GeV cm -3 . This goal has 
already been achieved by the DAMA experiment p[, which has reported the indication of 
an annual modulation effect in its counting rate, compatible with values of the quantity 

n ^ (nuclcon) • , i 

Scalar 111 the ran g e 

3 • 1(T 9 nbarn < fl?*^ < 1 • 1CT 8 nbarn, (1) 



for values of the WIMP mass, correlated with p x 3 cr^i^ on \ which extend over the range 30 

GeV < m x < 130 GeV 0. The region in the m x - P x ' 3 ^sc^ax° plane, singled out by the 
DAMA experiment at 2-cr C.L., is the one depicted in Fig. 2 of Ref. || and therein (and 
here) denoted by R m . Let us also notice that, taking into account the uncertainties in the 
local total dark matter density: 0.1 GeV cm -3 < pi < 0.7 GeV cm" 3 ]7j], when one assumes 
that a single WIMP candidate saturates pi, the range of Eq.(^) implies for cTg™ 1 ^ I r eon ' ) : 

1 ■ 10- 9 nbarn < < 3 ■ 10~ 8 nbarn, (2) 



Another experiment of WIMP direct detection which is now entering the upper part of the 
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range of Eq.(|T]) for Py'V^^ * s ^ ne CDMS experiment 

Once a given range for Px^icailr^ * s singled out by an experiment, what are the im- 
plications for specific particle candidates? We addressed this question in Refs. |9||| in the 
case of the DAMA data, under the assumption that the reported indication of an annual 
modulation is interpreted in terms of an effect due to relic neutralinos. We derived a num- 
ber of features for the neutralino cosmological properties, by selecting the supersymmetric 
configurations on the basis of the range of Eq. ([[]) , appropriately correlated with the range of 
m x . In particular, we proved that this range for p^' 3 ^™^ 01 ^ is compatible with a neutralino 
as a sizeable component of dark matter. In the derivation of this property, the role played by 
the correlation between cjg™^ 011 ^ and the neutralino-neutralino annihilation cross-section, 
<j ann , is crucial. In fact, these two cross-sections are normally both increasing or decreasing 
functions of the supersymmetric parameters. Since the neutralino relic abundance Q x h 2 (fl x 
is the neutralino cosmological density and h is the Hubble parameter in units of 100 km s _1 
Mpc -1 ) is roughly inversely proportional to a ann lower bounds on (Xg™]^ entail upper 
bounds for the relic abundance. It is remarkable that, as mentioned above, the range of 
Eq.(l), singled out by the indication of annual modulation, is compatible with a neutralino 
relic abundance of cosmological interest |6]||. 

From the previous discussion it is clear that one of the crucial ingredients in the derivation 
of the neutralino relic abundance from the results of direct detection experiments is ff^ 6011 ' . 
But, actually, how accurately can this quantity be evaluated at present? The point is that 
^scalar* usually takes dominant contributions from interaction processes, where neutralinos 
and quarks (inside the nucleon) interact by exchange of Higgs particles or squarks, and 
the relevant couplings are still plagued by sizeable uncertainties, related to hadron physics, 
which have not yet received satisfactory answers. Indeed, the Higgs-quark-quark or squark- 
quark-neutralino couplings involve the use of quark masses and quark scalar densities inside 
the nucleon, i.e. quantities which are subject to large uncertainties. 



We pointed out this problem in Ref. ||10|| , where we showed that, by taking two different 
determinations p"T| , |T2f of the pion-nucleon sigma term, <j w n, and of the fractional strange 
content in the nucleon y (see later on for the definitions of these two quantities), these 
uncertainties may affect cr^^° n ^ approximately by a factor of 3. This point, subsequently 



also recognized in Ref. JIB], was usually overlooked in subsequent papers, though the afore 



mentioned uncertainties still persist. The common practice prevailed of employing, for the 
relevant couplings, some standard values, which became consolidated more because of their 
reiterated use, rather than because of confirmation by more refined theoretical analyses. 

Actually, the quantities a^N and y have recently been the object of various other calcu- 
lations, based mainly on chiral perturbation theory and on QCD simulations on a lattice; 
however, the present situation, which is schematically revised in the next section, is still 
far from being clear. Thus, the Higgs-quark-quark and squark-quark-neutralino couplings 
are still plagued by significant uncertainties, of the order of those pointed out in Ref. |10 



In the present paper we wish to take up this problem again and, moreover, to address the 
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question about how sensitive are the neutralino cosmological and astrophysical properties, 
as derived from direct detection measurements, to the uncertainties which still affect the 
neutralino-nucleon cross-section. 

The theoretical framework employed here is the Minimal Supersymmetric extension of 
the Standard Model (MSSM), with the specifications given in Refs. to which we refer 
for all relevant details. Our results will mainly be given in terms of scatter plots, obtained 
by scanning the supersymmetric parameter space over the grid defined in Refs. Of 
particular importance for the properties discussed in the next section is the role of the 
neutral Higgs bosons: the two CP-even ones, h and H, and the CP-odd one, A. h and H 
are the main mediators of the coherent neutralino-nucleus cross-section, A is important in 
the neutralino-neutralino annihilation channels. The neutralino is defined as the lowest- 
mass linear superposition of photino (7), zino (Z) and the two higgsino states {H°, 

X = ai7 + a 2 Z + a 3 if ° + a^H^ . (3) 



To classify the nature of the neutralino, we define a parameter P = af + af; hereafter 
the neutralino is called a gaugino, when P > 0.9, is called mixed when 0.1 < P < 0.9 and a 
higgsino when P < 0.1. 



II. NEUTRALINO NUCLEUS ELASTIC CROSS SECTION 



We recall that the neutralino-nucleon scalar cross-section is given by 



cr, 



(nucleon) 
scalar 



7T 



mi 



m 



(4) 



where the two first terms inside the brackets refer to the diagrams with h- and if -exchanges 
in the t-channel (the A-exchange diagram is strongly kinematically suppressed and then 
omitted here) jnj and the third term refers to the graphs with squark-exchanges in the 
s- and u-channels [0. The mass m re d is the neutralino-nucleon reduced mass. Here, for 
simplicity, we explicitly discuss only the Higgs-mediated terms, then we do not report the 
expressions for the squark propagator Pg. and for the couplings Ag., Bg., which may be 
found in Ref. ||. However, arguments analogous to the ones given below hold also for the 
squark-exchange terms, which are actually included in the numerical results reported in this 
paper. 



The quantities F h)H and Ih,H are defined as follows 

F h = (— a,\ sin W + a 2 cos 9w) sin a + a 4 cos a) 
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Fh = (— ai sin 6w + «2 cos ^vk) (03 cos a — CI4 sin a) 
4^ = £^m g (iV|gg|iV>. (5) 



The matrix elements < N\qq\N > are meant over the nucleonic state. The angle a rotates 
h[ ^ and Hq^ into h and if, and the coefficients are given by 

k t-t yP e = cos a/ sin /?, fcf_ type = - sin a/ sin (3 

kd-type = ~ Sil1 «/ COS ft ^f-type = ~ C °S «/ COS /? (6) 

for the up-type and down-type quarks, respectively 

The quantities m q (N\qq\N)'s for the light quarks u,d,s may conveniently be expressed 
in terms of the pion-nucleon sigma term 

o'nN = 7^ ( m « + m d) < N\uu + dd\N >, (7) 

the fractional strange-quark content of the nucleon 

< iV|ss|iV > 

and the ratio r = 2m s /(m u + mj. In fact, assuming isospin invariance for quarks u and cf, 
one has 

m u < N\uu\N > ~ m d < N\dd\N >~ ^a nN (9) 
m s < N\ss\N > ~ -rye^N- (10) 

For the heavy quarks c, b, t, using the heavy quark expansion jlBf , one derives 

m c < iV|cc|iV > ~ m b < iV|6&|iV >~ m t < N\it\N >~ 
2 

~ 27 



1 , 

m N - (1 + -ryjo-irN 



(11) 



where is the nucleon mass. The quantities i^n can then be rewritten as 

-typedu + ™d—type9d> 

where 



f/i,H — k u '_ type g u + k d !_ t gj, (12) 



4 . 19 1 , 2 . 23 25 

g u = —{m N + ~ -rya^N), gd = ^{m N + —a nN + —rya nN ). (13) 

We turn now to the values to be associated to <J n N, y and r. 
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A. Pion nucleon sigma term. 



The quantity a^N ma Y be deduced phenomenologically from measurements of the pion- 
nucleon scattering; however, its derivation from the experimental data is rather involved. 
The customary procedure is to go through the following steps (see, for instance, Ref. |17||): 

i) By use of phase-shift analysis and dispersion relations, from the experimental data of 
low-energy pion-nucleon scattering, one derives the quantity 

S CD = £(* = 2ml) = JXn(s = m 2 N} t = 2ml), (14) 

where s and t are standard Mandelstam variables, m v is the pion mass, /„. is the pion-decay 
constant and T nN is the (Born-subtracted) pion-nucleon scattering amplitude, calculated 
at the so-called Cheng-D ashen point. 

ii) Modulo terms of order < 1 MeV, which may be safely neglected, one has 

Z CD ~a 7TN (t = 2ml) J (15) 

where a^^it) is the nucleon scalar form factor, defined as 

crMt = (p' - P?) =< N(p')\ mu+ 2 md (uu + dd)\N(p) > . (16) 



iii) The evolution of G n N(t), as a function of the momentum transfer from t = 2ml to 
- 0, provides the value of g^n = oVtv^ = 0) : 

a nN = a nN (t = 2m 2 n ) - Ag. ( 17 ) 



Now, it turns out that the determinations of Tjcd and Ag suffer from sizeable uncertain- 



ties. 



For Ticdi apart from older calculations, we have, for instance: T,cd =56-r 72 MeV |T8 



and S CD =59^- 62 MeV JT^JT^]. A value of S CD on the high side (~ 72 MeV) is also favoured 
by a more recent evaluation pTj. Thus, one could tentatively consider the range 



E CD = 56 -f- 72 MeV. 



(18) 



Calculations of Ag by dispersion relation techniques provide [|19[ (see also [plf] ): 
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Act=15.2±0.4 MeV. 



(19) 



These values are much larger than the ones obtained with chiral perturbation theory at 
leading order Act ~ 7.5 MeV [E3j (see Ref. [33] for a possible explanation of this discrepancy). 
Furthermore, it has to be noted that Act as deduced from lattice calculations J3JJ is given 
by Acx=6.6±0.6 MeV, thus by values sizably smaller than those of Eq.flTDl). 



Combining the range for Sod i n Eq.(|l~8"D and the value AS=15 MeV, we obtain for ct^n 
the range 41 MeV< a n N <57 MeV. This may be compared to the results derived in Ref. 
pl| using heavy quark chiral perturbation theory: 38 MeV< ct^n <58 MeV and in Ref. |26| 
by lattice calculations: 40 MeV< a^N <60 MeV. Thus we conclude that the value of a nN 
is still considerably uncertain; the previous results only indicate some convergence towards 
the range 40 MeV < a n N < 60 MeV, with an upper extreme which might be even higher 
(~ 65 MeV), should one take the chiral perturbation theory result (Act ~ 7 MeV), instead 
of the dispersion-relation result (Eq.(|T9"|)). Finally, we wish to notice that recent results 
from higher order chiral perturbation calculations provide a large value for ct^n'- a n ]y = 70 
MeV [27]; furthermore, use of a new pion-nucleon phase-shift analysis P5 |, instead of the 
standard one [^] would make cr^N to rocket to a value larger than 200 MeV 



B. Strange quark content of the nucleon. 



A standard way to evaluate the quantity y defined in Eq.([|) is to express it in terms of 
a nN and of the quantity cr defined as 



ct = -(m u + m d ) < N\uu + dd — 2ss|iV >, 
2 



(20) 



i.e. 



go 



(21) 



Actually, cr is a quantity related to the size of the SU(3) symmetry breaking and, as such, 



may be calculated either from the octect baryon masses: o"o —33 MeV |TTJ or with chiral 



perturbation theory: cr =35±5 MeV |22j, cr =36±7 MeV [31]. For definiteness, we take 



a = 30 40MeV. 



(22) 



Thus, we have, for instance: 0< y < 0.25 for a 7r Ar=40 MeV, 0.11< y < 0.33 for <r nN =4:5 
MeV, 0.33< y < 0.50 for a nN =Q0 MeV, and 0.38< y < 0.54 for a wN =65 MeV. 



7 



We wish to remark that, by general physical arguments, one would expect for y a some- 
what small value, i.e. y ~ 0.2-0.3; however, apart from the results of the previous deriva- 
tions, which allow values of y up to y ~ 0.5, also lattice results seem to favour large values: 
y=0.36±0.03 |24|, y=0.66±0.15 ]26| (even reducing this latter value by ~ 35% as suggested 



in p4"|, one would still have y ~0.4-0.5). 



C. Mass ratio r = 2m s /(m u + m^). 



As we have seen at the beginning of this section, besides a w N and y a third ingredient is 
necessary for the evaluation of m q (N\qq\N) in the case of the strange quark and the heavy 
ones: the mass ratio r = 2m s / (m u + rud)- 

The standard derivation of this ratio is based on chiral perturbation theory. Lowest 
order formulae (corrected for electromagnetic effects) which give the mass ratios m u /m^ 
and m s / rrid m terms of the physical masses of the K mesons, entail, for m s j (m u + rrid), the 



canonical value 32 



26. (23) 



The inclusion of next-to-leading order contributions in the chiral expansion lead to the 
determination [ |33| 1 



24.4 ± 1.5. 



(24) 



Use of mass ratios in the evaluation of r makes these results independent of the renormaliza- 
tion scale. However, the validity of the chiral perturbation method relies on the hypothesis 
that the quark condensate is the leading order parameter of the spontaneously broken sym- 
metry. 

Other methods, most notably QCD sum rules and lattice simulations of QCD, are capable 
of providing evaluations of individual quark masses (not only of their ratios), though these 
derivations still suffer from large uncertainties. For the u and d quark masses we can quote 
the following results derived from QCD sum rules: 

(m u + m d )(l GeV) = 12.0 ± 2.5 [34] (25) 
(m u + m d )(l GeV) = 15.5 ±2.0 ||. (26) 



Even more spread are the values deduced with this method by various authors for the 
s-quark mass [p6|-|40[|. A combination of all these derivations: m s (l GeV) = 170 ± 50 MeV 
PH evidentiates how large is the uncertainty in the estimate of m s . For sake of comparison, 
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we can quote some values for m s as derived from lattice QCD: m s (l GeV) = 155 ± 15 MeV 
jg-H and from the r hadronic width: m s (l GeV) = 193 ± 59 MeV g| and m s (l GeV) = 



200 ± 70 MeV M 



If, following Ref. [|3|], one uses 

m s (l GeV) = 175 ± 25 MeV 



(27) 



(which combines the results of Refs. fl3"6|j3~?|j ) and takes m u + rrid = 12.0 ± 2.5 (see Eq.(|25] 
one obtains 



r = 29 ± 7. 



(28) 



This result may be considered as representative of the uncertainty currently affecting the 
mass ratio r. 



D. Size of the Higgs quark couplings. 



From Eqs.(|9l, |i0yTTD and the ranges for a n N, V and r discussed in Sections [II A| - |II C| , we see 
that the quantities m q (N\qq\ N)'s are indeed affected by large uncertainties. In particular, 
the quantity m s (N\ss\N) , is uncertain roughly by a factor of 3-4. This has quite significant 
consequences on the size of CTg™^ 0011 ^ , because m s < N\ss\N > is the most important term 
among the m q < N\qq\N >'s K7L unless tan/? is very small (see Eqs. (|6|),(|T2]),(p!3|)). 



The values of the m q {N\qq\N) , s for a few sets of values for a^N, V and r are given in 
Table I. The first three sets are representative of values currently employed in the literature. 
In particular the third set (denoted here as set 1) is the one used in our previous papers 
IIU on the analysis of the DAMA annual modulation data. 

Set 2 and set 3 are representative sets of values which are employed here to illustrate: i) to 
which extent the size of (r^ilr"^ may be increased, within the afore mentioned uncertainties, 
and ii) which are the ensuing implications for the neutralino cosmological properties, when 
these are derived from experimental data of WIMP direct detection. Thus, sets 2 and 3 
are meant only to be two possible representative set of values, chosen within the current 
uncertainties and capable of providing large values of the quantity m s < N\ss\ N >. Because 
of the correlations, previously discussed, among o"g°^ 1 r con ' ) , a ann and Vt x h 2 , sets 2 and 3 are 
expected to provide sizeable values for the neutralino relic abundance. 

A more systematic analysis of the implications for relic neutralinos of uncertainties in 
the Higgs-quark couplings will be feasible, only when a consistent field of variation for 
the correlated quantities ovjv, y and r will emerge from a thorough and coherent QCD 
investigation of the problem. 
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III. RESULTS AND CONCLUSIONS 



Let us now turn to the presentation of our results. In Fig.l.a (l.b) we give the ratio of 
the cross-section cr^i^° calculated with set 2 (set 3), to cT^ alar eon ' ) calculated with set 1. 
We see that for (o"scaiar 60n ' ) )set i in the range of Eq. (0) most configurations cluster around the 

i / (nucleon) \ ; / (nucleon) \ / (nuclcon) \ / / (nucleon) \ r mi 1 

values « alar ; )set2/« a i ar )set i ^ 3, (a s ^ calar ; )set3/« a iar )set i ^ 5. Thus, we have a 
sizeable increase in the cross-section, when sets 2 or 3 are used instead of our set of reference 
(set 1), which was utilized in Refs. ||,|6|]. 

Fig. 2 displays the plot of cr s ^ alar eon ' ) versus Q x h 2 , in case of set 1 (Fig. 2a) and of set 2 
(Fig. 2b). The two horizontal dashed lines delimit the range of the neutralino- nucleon cross 
section defined in Eq. (|2|) . The solid vertical lines delimit the cosmologically interesting range 
0.01 < Q x h 2 < 0.7. The two vertical dashed lines delimit the range: 0.02 < £l x h 2 < 0.2, 
which represents a particularly appealing interval, according to the most recent observations 



and analyses |48|Ji9| . The overall shape of these scatter plots reflects the anticorrelation 
between crg^ l ala 1 r eon ' ) and Q x h 2 . The most relevant feature of these plots is the boundary on the 
top-right side. It provides the maximal values allowed for the relic density, for neutralino- 
nucleon cross-sections in the range of Eq. (^). Notice how the boundary extends into the 
region of cosmo logical interest more markedly in case of set 2 than in case of set 1 . This is a 
first manifestation of the fact that the DAMA annual modulation effect is compatible with 
a neutralino of cosmo logical interest, to an extent which is even larger than the one singled 
out in our previous papers where only the representative set 1 was employed. (Notice 
that in the present analysis, although the overall range of variation of the susy parameters is 
the same as in Ref . |9],|6| , we have optimized the numerical scanning of the parameter space 
in order to increase the density of configurations which fall in the region of main interest.) 

These cosmological properties are further displayed in Fig. 3, which depicts the scatter 
plots of p x versus Q x h 2 . Here the two horizontal lines delimit the physical region for the 
total local density of non-baryonic dark matter: 0.1 GeV cnr 3 < p x < 0.7 GeV cm -3 ; the 
two slant dot-dashed lines delimit the band, where linear rescaling procedure for the local 
density is usually applied At variance with the previous scatter plots of Figs. 1-2, 

which refer to a generic scanning of the supersymmetric parameter space, constrained only 
by accelerator bounds, as discussed in Refs. the scatter plots of Fig. 3 display the susy 
configurations singled out by the DAMA annual modulation data. These plots are obtained 
with the following procedure: 

. \ i , i r (nucleon) n / (nucleon) i r (nucleon) n n , ,1 , 

i) p x is evaluated as p x = [p x < alar ; ]i? m /^ calar ; , where [p x c^ calar >\ Rm denotes the set 
of experimental values of p x Scalar ° n ' > inside the DAMA annual modulation region R m and 
^scalar' 011 '' i s calculated with Eq.(|J). 

ii) To each value of p x , which then pertains to a specific supersymmetric configurations, 
one associates the corresponding value of Q x h 2 , calculated as in Ref. |50 . 
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Therefore, with this procedure, we determine the values of p x which, for each calculated 
CT scXx° > satisfy the DAMA annual modulation data. 



Fig. 3 shows that the set of supersymmetric configurations selected by the DAMA data 
has a significant overlap with the region of main cosmological interest: Q x h 2 > 0.02 and 0.1 
GeV cm" 3 < p x < 0.7 GeV cm -3 . The extent of this overlap is increasingly larger for set 
2 and set 3 of Table I. By way of example, for set 3 one has that, at p x = 0.3 GeV cm -3 , 
Q x h 2 may reach the value 0.3. Therefore, these results reinforce our conclusions of Ref. ||||, 
i.e. that the DAMA annual modulation data are compatible with a neutralino as a major 
component of dark matter, on the average in the Universe and in our Galaxy. 

The same figure shows that different situations are also possible. Specifically, for config- 
urations which fall inside the band delimited by the slant dot-dashed lines, the neutralino 
would provide only a fraction of the cold dark matter both at the level of local density and at 
the level of the average Q, a situation which would be possible, for instance, if the neutralino 
is not the unique cold dark matter particle component. On the other hand, configurations 
above the upper dot-dashed line and below the upper horizontal solid line would imply a 
stronger clustering of neutralinos in our halo as compared to their average distribution in 
the Universe. This situation may be considered unlikely, since in this case neutralinos could 
fulfill the experimental range for p x , but they would contribute only a small fraction to the 
cosmological cold dark matter content. Finally, configurations above the upper horizontal 
line are incompatible with the upper limit on the local density of dark matter in our Galaxy 
and must be disregarded. 



Finally, we wish to add the following comments: 



i) To show the effect of the current uncertainties in the Higgs-quark couplings, a few 
representative sets of values for the relevant quantities (a^N, V and r) were employed, as 
illustration of the implications for relic neutralinos. We have explicitly discussed situations 
which allow sizeable values of the neutralino relic abundance. It is clear that, within the same 
uncertainties, other sets of values of the relevant parameters would entail relic neutralinos of 
much less appealing cosmological interest. A more precise evaluation of the effects implied 
by the mentioned uncertainties will require a significant breakthrough in the understanding 
of the QCD and of other hadron aspects involved in the problem. What is lacking at present 
is consistency among the various determinations in a n N, V and r. 

ii) The relevance of the implication of the mentioned uncertainties in converting infor- 
mation from measurements of neutralino detection to neutralino cosmological properties, 
discussed here for experiments of direct detection, applies also to indirect experiments at 
neutrino telescopes, when one looks at upgoing muons from the center of the Earth. In fact, 
in this case the size of the expected signal depends on the neutralino capture rates by the 
Earth, and then mainly on p^^^^' 

iii) A number of experimental means, meant to search for the presence of relic neutralinos 
in our galaxy, look for signals (antiprotons, antideuterons, positrons, diffuse gamma-rays, 
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gamma lines) [|TJ originated by neutralino-neutralino annihilation in the halo. Thus, the 
relevant cross-section involved in this case is cr ann . Therefore, a word of caution is in order 
here, for the case when the expectations for these indirect signals are based on supersym- 
metric configurations derived from direct detection data. In the estimate of these indirect 
signals, one should include the uncertainties discussed in the present paper. 
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TABLES 



(MeV) 


y 


r 


m qi < N\qiqi\N > 
(MeV) 


m s < N\ss\N > 
(MeV) 


m h < N\hh\N > 
(MeV) 


Ref. 


45 


0.28 


25 


23 


158 


55 


|1C] 


45 






27 


131 


56 


011 


45 


0.33 


29 


23 


215 


50 


set 1 


60 


0.50 


29 


30 


435 


33 


set 2 


65 


0.50 


36 


33 


585 


21 


set 3 



TABLE I. Values of the matrix elements (N\qq\N) of the quark scalar densities in the nucleon 
times the quark masses m q , for a few sets of values of the pion-nucleon sigma term ct^n, the 
fractional strange-quark content of the nucleon y and the quark mass ratio r = 2m s /(m u + m^). 
qi stands for light quarks, s is the strange quark and h = c,b,t denotes heavy quarks. For the light 
quarks, we have defined m qi < N\q~iqi\N > = \\m u < N\uu\N > +md < N\dd\N >]. 
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FIGURE CAPTIONS 



FIG. la. Ratio of the neutralino-nucleon scalar cross-section (^catex ) > calculated 
with the parameters of set 2 defined in Table I, to (cr^aiaT ) _ t : calculated with set 1, as a 
function of (a. 



(nucleon) \ 
scalar J 



set 1 



FIG. lb. The same as in Fig. la, with set 3 instead of set 2. 



FIG. 2a. Neutralino-nucleon scalar cross-section cr s l °^ r eon , calculated using the param- 
eters of set 1, as a function of the neutralino relic abundance Q x h 2 . The two horizontal 
dashed lines delimit the range of the neutralino-nucleon cross section defined in Eq.(|2|). The 
two solid vertical lines delimit the interval of cosmological interest. The two vertical dashed 
lines delimit the preferred band for cold dark matter. The shaded region is cosmologically 
excluded on the basis of present limits on the age of the Universe. 



FIG. 2b. The same as in Fig.2a, for the parameters of set 2. 



FIG. 3a. Neutralino local density p x derived by requiring that p x a S caiax ^ aus inside the 
experimental DAMA region R m , plotted against the neutralino relic abundance Q x h 2 . The 
quantity <J^^ on ^ is calculated with the parameters of set 1. The two horizontal lines delimit 
the physical range for the local density of non-baryonic dark matter. The two solid vertical 
lines delimit the interval of Q x h 2 of cosmological interest. The two vertical dashed lines 
delimit the preferred band for cold dark matter. The two slant dot-dashed lines delimit the 
band where linear rescaling procedure is usually applied. The shaded region is cosmologically 
excluded on the basis of present limits on the age of the Universe. Different symbols identify 
different neutralino compositions: circles stand for a higgsino, crosses for a gaugino and dots 
for a mixed neutralino. 



FIG. 3b. The same as in Fig.3a, for the parameters of set 2. 



FIG. 3c. The same as in Fig. 3a, for the parameters of set 3. 



16 



FIGURES 




FIG. la - Ratio of the neutralino-nucleon scalar cross-section (cicaiar° n ^) > calculated 
with the parameters of set 2 defined in Table I, to (cr^^° n ^ t i calculated with set 1, as a 

r ,. r ( (nucleoli) \ 

function of (^ calar 'J^. 
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FIG. lb - The same as in Fig. la, with set 3 instead of set 2. 



18 



A, Esj-ttlnci, r. Donate, tf, rgmenjci. B, Baopel 




FIG. 2a - Neutralino-nucleon scalar cross-section cr s °^ r eon , calculated using the param- 
eters of set 1, as a function of the neutralino relic abundance Q x h 2 . The two horizontal 
dashed lines delimit the range of the neutralino-nucleon cross section defined in Eq. (|2|) . The 
two solid vertical lines delimit the interval of cosmological interest. The two vertical dashed 
lines delimit the preferred band for cold dark matter. The shaded region is cosmologically 
excluded on the basis of present limits on the age of the Universe. 
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FIG. 3a - Neutralino local density p x derived by requiring that Px^scaiar f a ^ s inside the 
experimental DAM A region R m , plotted against the neutralino relic abundance Vl x h 2 . The 
quantity CscaiaT™^ is calculated with the parameters of set 1. The two horizontal lines delimit 
the physical range for the local density of non-baryonic dark matter. The two solid vertical 
lines delimit the interval of Q x h 2 of cosmological interest. The two vertical dashed lines 
delimit the preferred band for cold dark matter. The two slant dot-dashed lines delimit the 
band where linear rescaling procedure is usually applied. The shaded region is cosmologically 
excluded on the basis of present limits on the age of the Universe. Different symbols identify 
different neutralino compositions: circles stand for a higgsino, crosses for a gaugino and dots 
for a mixed neutralino. 
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FIG. 3b - The same as in Fig. 3a, for the parameters of set 2. 
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FIG.3c - The same as in Fig. 3a, for the parameters of set 3. 
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